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Abstract 

In a genome-wide association study, we have previously identified and per- 
formed the initial replication of three novel susceptibility loci for cervical can- 
cer: rs9272143 upstream of HLA-DRBl, rs2516448 adjacent to MHC class I 
polypeptide-related sequence A gene (MICA), and rs3117027 at HLA-DPB2. 
The risk allele T of rs25 16448 is in perfect linkage disequilibrium with a frame- 
shift mutation (A5.1) in MICA exon 5, which results in a truncated protein. To 
validate these associations in an independent study and extend our prior work 
to MICA exon 5, we genotyped the single-nucleotide polymorphisms at 
rs9272143, rs25 16448, rs3117027 and the MICA exon 5 microsatellite in a 
nested case-control study of 961 cervical cancer patients (827 carcinoma in situ 
and 134 invasive carcinoma) and 1725 controls from northern Sweden. The C 
allele of rs9272143 conferred protection against cervical cancer (odds ratio 
[OR] = 0.73, 95% confidence interval [CI] = 0.65-0.82; P = 1.6 x 10"^), 
which is associated with higher expression level of HLA-DRBl, whereas the T 
allele of rs2516448 increased the susceptibility to cervical cancer (OR = 1.33, 
95% CI = 1.19-1.49; P = 5.8 x 10"^), with the same association shown with 
MICA-A5.1. The direction and the magnitude of these associations were consis- 
tent with our previous findings. We also identified protective effects of the 
MICA-A4 (OR = 0.80, 95% CI = 0.68-0.94; P = 6.7 x 10"^) and MICA-A5 
(OR = 0.60, 95% CI = 0.50-0.72; P = 3.0 x lO"**) alleles. The associations with 
these variants are unlikely to be driven by the nearby human leukocyte antigen 
(HLA) alleles. No association was observed between rs3117027 and risk of 
cervical cancer. Our results support the role of HLA-DRBl and MICA in the 
pathogenesis of cervical cancer. 
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Introduction 

Worldwide, cervical cancer is the third most common 
cancer and the second most frequent cause of cancer 
deaths among women, and resulted in an estimated 
530,000 new cancer cases and 275,000 deaths in 2008 [1]. 



In many low-income countries, cervical cancer is the most 
common cancer and the leading cause of cancer-related 
death among women [ 1 ] . Cervical cancer and its precursor 
lesions, cervical intraepithelial neoplasia (CIN) are caused 
by persistent infection with high-risk human papillomavi- 
rus (HPV), where CIN III is considered the same as carci- 
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noma in situ (CIS) [2]. During their lifetime, many 
women will become infected with HPV, but only a minor- 
ity will develop CIN or cervical cancer. Consequently, 
other factors, for example, host genetic factors, play an 
important role in both the persistence of infection and 
progression to cancer [3, 4]. 

We have recently performed the first genome-wide asso- 
ciation study (GWAS) of cervical cancer and identified 
three independent novel loci within the major histocom- 
patibility complex (MHC) region at 6p21.3 that influence 
susceptibility to cervical cancer in a Swedish population. 
The first is located between HLA-DRBl and HLA-DQAl 
(rs9272143; odds ratio [OR] = 0.67, 95% confidence inter- 
val [CI] = 0.62-0.72 for C allele; P = 9.3x10"^"); the sec- 
ond is adjacent to the MHC class I polypeptide-related 
sequence A gene (MICA) (rs2516448; OR = 1.42, 95% 
CI = 1.31-1.54 for T allele; P = 1.6 x 10""*); and the 
third at HLA-DPB2 (rs3117027; OR = 1.25, 95% 
CI = 1.15-1.35 for A allele; P = 4.9x10"**) [5]. The associ- 
ations observed for these three new loci were found to be 
statistically independent of previously known associations 
with the human leukocyte antigen (HLA) alleles/haplotypes 
[5]. The transmembrane domain (TMD) of MICA encoded 
by exon 5 harbors a variable number of GCT repeats, which 
encode 4, 5, 6, or 9 alanine (Ala) residues (alleles designated 
A4, A5, A6 or A9, respectively). Additionally, the A5.1 allele 
(rs67841474) contains an extra guanine (G) insertion after 
two GCT triplets, which causes a frameshift mutation 
resulting in a premature stop codon that, in turn, truncates 
10 amino acids of the TMD as well as the hydrophobic 
cytoplasmic tail [6]. The risk allele T of rs2516448 was 
found to be in perfect linkage disequilibrium (LD) (D' = 1, 

= 1) with A5.1 and cervical neoplasia patients carrying 
the A5.1 allele have less membrane-bound MICA in their 
lesions [5]. 

In our initial GWAS, we were able to replicate the 
effect of the three susceptibility loci in a second cohort 
from southern and middle Sweden. However, validation 
of GWAS findings in multiple cohorts is necessary in 
order to report genotype-phenotype associations. The 
new susceptibility loci for cervical cancer require further 
investigation in a large sample size. It is also important 
to extend our prior work to MICA exon 5 microsatel- 
lite polymorphism and evaluate effect modification by 
age of onset and tumor stage. Therefore, we investi- 
gated the association between single-nucleotide poly- 
morphisms (SNPs) of rs9272143, rs2516448, and 
rs3117027 as well as MICA exon 5 microsatellite poly- 
morphism and risk of cervical cancer, in a large nested 
case-control study of 961 incident cervical cancer 
patients (827 CIS and 134 invasive carcinoma) and 
1725 cancer- free controls from the Vasterbotten County 
in northern Sweden. 



Extension of Findings in Cervical Cancer GWAS 

Material and Methods 
Study population 

Eligible women for the study were defined as Vasterbotten 
County resident in northern Sweden who had at least one 
cytologicaUy normal cervical smear and who had no prior 
operative treatment of the cervix. Linkage between the 
cytology registry and the Swedish Cancer Registry from 
1961 identified 832 patients with CIS and 134 patients 
with invasive cervical cancer diagnosed after the sampling 
date of a normal smear. Controls were women in the 
study base who did not develop cervical cancer before the 
time point of diagnosis of a corresponding case. For each 
CIS case, two population-based controls were selected, 
matched for age of subject (±5 years) when the sample 
was collected. For each invasive cervical cancer case, one 
population-based control was selected, matched for age of 
subject (±5 years) when the sample was collected. A writ- 
ten informed consent was obtained from each participant 
and this study was approved by the Institutional Review 
Board (IRB) of the Umea University. Genomic DNA was 
extracted from the buffy coat using standard phenol- 
chloroform extraction protocol. In total, DNA samples 
from 827 women with CIS and 1591 matched healthy 
controls, and 134 women with invasive cervical cancer 
and 134 matched healthy controls qualified for genotyp- 
ing. The study population was not included in the previ- 
ous GWAS study. 

Genotyping assay 

Single-nucleotide polymorphisms of rs9272143 and 
rs3 117027 were genotyped with template-directed dye- 
terminator incorporation with fluorescence polarization 
detection (FP-TDI) (Tecan, Mannedorf, Switzerland) and 
rs25 16448 was genotyped using the TaqMan assay (Applied 
Biosystems, Foster City, CA) . The information of the prim- 
ers and probes is described in Supplementary Table SI. The 
polymerase chain reaction (PCR) amplification of the 
MICA microsatellite alleles of exon 5 was carried out using 
a 5' end fluorescently (6-FAM)-labeled reverse primer and 
a forward unlabeled primer. The primer sequences for the 
MICA microsatellite were previously reported [7[ and are 
described in Supplementary Table SI. The PCR products 
were mixed with Hi-Di Formamide and GeneScan 500 
ROX size standard and separated on an ABI 3730x1 DNA 
Analyzer (Applied Biosystems). Different alleles were anno- 
tated using GeneMapper 4. 1 software (Applied Biosystems, 
Foster City, USA) based on the size of the PCR products. 
Eight percent of the samples were selected for repeat geno- 
typing as duplicates, yielding a reproducibility of 100%. 
Genotype success rate was >98.21%. 
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Table 1. Selected demographic clnaracterlstics of study subjects. 



Cases Controls 
A/2 (%) A/2 (%) 



Age (years)' 






<36 


478 (49.33) 


851 (49.74) 0.22 


>36 


483 (50.67) 


874 (50.26) 


Mean± SD 


36.86 ± 9.05 


36.43 ± 8.86 


Study design 






Study 1 


827* (86.06) 


1591 (92.23) 


Study 2 


134= (13.94) 


134 (7.77) 


Total 


961 


1725 


SD, Standard error. 






'The median age 


is 36 years for both 


cervical cancer patients and 



control subjects. 
^Number of samples. 

^Difference in age between cervical cancer patients and control sub- 
jects was evaluated by using f-test. /'value is two-sided. 
"Number of subjects with carcinoma in situ. 
^Number of subjects with invasive carcinoma. 

Statistical analyses 

Differences in age between the cervical cancer cases and 
controls were evaluated by using t-test. Goodness-of-fit to 
the Hardy-Weinberg equilibrium (HWE) expectation in 
control subjects was assessed by a chi-square (x^) test 
implemented in PyPop for each locus [8]. The association 
between each genetic variant and disease risk was estimated 
by the odds ratio (OR) and 95% confidence interval (CI) 
per allele, assuming a log-additive genetic model using 
unconditional logistic regression adjusting for age at 
recruitment and study design (CIS vs. invasive carcinoma) 
in 961 cervical cancer patients and 1725 controls. 

For genetic variants that showed evidence of association, 
the heterozygous and homozygous carriers of the variant 
allele were compared with the noncarriers, respectively. 
Analyses were also conducted after stratifying for study 
design or tumor stage (CIS vs. invasive carcinoma) and age 
at recruitment (age<36 and age>36), adjusting for age at 
recruitment and study design (CIS vs. invasive carcinoma) 
when necessary. Heterogeneity of odds ratios across the 
stratification groups was assessed using the Cochran Q test. 
The potential interaction between each genetic variant and 
study design (CIS vs. invasive carcinoma) was evaluated 
with the interaction model by additionally including inter- 
action term between genotypes and study design. Statistical 
analyses were all performed using SAS 9.3 software (SAS 
Institute, Cary, NC), with two-sided tests. A Bonferroni 
correction for multiple tests was applied and gave a P value 
of 7.1 X 10^^ as the cutoff for statistical significance based 
on seven independent genetic variants tested (SNP 
rs2516448 is in perfect LD with A5.1 of MICA, hence they 
represent one independent test). 
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Results 

The characteristics of the cervical cancer patients and 
cancer-free controls enrolled in the study are described in 
Table 1. Overall, there was no significant difference in 
age between the cervical cancer patients and the control 
subject [P = 0.22), suggesting that matching based on age 
was adequate. The median age was 36 years for both 
cases and controls. Of the 961 cervical cancer patients, 
827 (86.06%) had a diagnosis of CIS and 134 (13.94%) 
of invasive carcinoma. 

Table 2 summarizes the estimates of the main effects 
for each SNP. Genotype frequency distributions in the 
control subjects were consistent with those expected from 
the HWE model for all SNPs (all P > 0.05). The variant 
allele T of rs25 16448 in the MHC class I region was sig- 
nificantly associated with increased risk of cervical cancer 
(OR = 1.33, 95% CI = 1.19-1.49; P = 5.8 x 10"^), 
whereas the variant allele C of rs9272143 in the MHC 
class II region was strongly associated with decreased risk 
of cervical cancer (OR = 0.73, 95% CI = 0.65-0.82; 
P= 1.6x10^^). Both the direction and magnitude of 
these associations were in accordance with our previous 
findings [5]. In contrast, there was no association between 
rs3117027 and risk of cervical cancer (OR = 0.99, 95% 
CI = 0.88-1.12 for variant allele A; P = 0.86). The LD 
between these three SNPs was very weak (r^ = 0) (Supple- 
mentary Table S2), consistent with previous study [5]. 

The allele frequencies of the MICA exon 5 microsatel- 
lite in cervical cancer patients and control subjects are 
shown in Table 3. Analysis showed that MICA-A5.1 (G 
insertion of rs67841474) had the highest frequency in 
both patients (60%) and control subjects (52%). In accor- 
dance with our previous finding [5], this microsatellite 
allele was in perfect LD (D' =1, = 1) with the risk 
allele T of rs2516448 in both cases and controls, and 
showed a comparable association with susceptibility to 
cervical cancer (OR = 1.34, 95% CI = 1.20-1.50; 
P = 3.8 X 10"^) as the rs2516448 T allele. In the overlap- 
ping 943 cervical cancer patients and 1683 cancer-free 
control subjects, the OR (95% CI) was 1.34 (1.20-1.51) 
for both A5.1 and the T allele of rs25 16448. In addition, 
significant protective effects were seen for MICA-A4 
(OR=0.80, 95% CI = 0.68-0.94; P = 6.7x10"^) and 
MICA-A5 (OR = 0.60, 95% CI = 0.50-0.72; 
P = 3.0 X 10^**). No correlation was found between 
rs9272143 and rs3117027 and alleles of MICA exon 5 
microsatellite (r^ = 0) (Supplementary Table S2). 

An allelic dosage effect on cervical cancer risk was 
observed for variant at rs9272143 in the MHC class II 
region and the MICA-A4, -A5 and -A5.1 alleles, when 
comparing the heterozygous and homozygous carriers of 
the variant allele with the noncarriers (Fig. 1). In particu- 
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Table 3. 


Association 


between MICA 


microsatellite and 


cervical 


can- 
















Allele count (allele 










frequency) 




Assnri^itinn'' 








Cases 


Controls 










(Total 


(Total 










alleles = 


alleles = 










1898)' 


3378)' 




p 




A4 


261 (0.14) 


566 (0.17) 


0.80 (0.68-0.94) 


6.7 X 


10-3 


AB 


179 (0.09) 


500 (0.15) 


0.60 (0.50-0.72) 


3.0 X 


10-*^ 


AB.1 


1138 (0.60) 


1771 (0.52) 


1.34 (1.20-1.50) 


3.8 X 


10"^ 


A5 


182 (0.10) 


290 (0.09) 


1.12 (0.92-1.36) 


0.25 




A9 


138 (0.07) 


251 (0.07) 


0.98 (0.79-1.22) 


0.87 





OR, odds ratio; CI, confidence interval. 

'Total allele counts in cervical cancer patients and control subjects, 
respectively. 

^Odds ratios and 95% confidence intervals for each allele in log-addi- 
tive model were derived from unconditional logistic regression adjust- 
ing for age at recruitment and study design (carcinoma in situ vs. 
invasive carcinoma). Two-sided P values correspond to the odds ratios. 

lar, MICA-A5.\ homozygotes had a 1.84-fold increased 
risk of developing cervical cancer (OR = 1.84, 95% 
CI = 1.46-2.32; P = 3.2 x 10"^), whereas UlCA-AS ho- 
mozygotes had nearly threefold protection against cervical 
cancer (OR = 0.34, 95% CI = 0.16-0.75; P = 7.0x10"'') 
as compared to noncarriers, respectively. No statistically 
significant heterogeneity was observed by age for any of 
the genetic variants. However, the variant of rs9272143 in 
the MHC class II region showed strong heterogeneity 
when stratifying by study design (or tumor stage) 
{P.het = 1 X 10"*), with little evidence for association in 
invasive carcinoma (OR = 1.38, 95% CI = 0.98-1.95; 
P = 0.07), although the number of subjects with invasive 
cancer was limited. By contrast, there was little evidence 
for heterogeneity by tumor stage for the MICA variants, 
and a statistically significant association with MICA-A5.1 
was observed in invasive carcinoma (OR = 1.46, 95% 
CI = 1.03-2.07; P = 0.03). Consistently, statistically signif- 
icant interaction was observed between tumor stage and 
rs9272143 (P = 4.4 x 10""), but not between tumor 
stage and MICA alleles (all P > 0.05). 

Discussion 

We replicated the associations of cervical cancer with 
rs9272143 located in the MHC class II region as well as 
with rs2516448 and MICA-A5.1 in the class I region iden- 
tified in our previous GWAS, with ORs of similar magni- 
tude to that previously reported [5]. We also identified 
protective effects of both the MICA-A4 and MICA-A5 
alleles against cervical cancer. None of these variants 
showed heterogeneity by age of onset. The association 
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rs9272143 

Overall 

TC versus TT 
CC versus TT 



Ca Co OR 



95%CI 



960 1725 0.73 0.65-0.82 
465 867 0.76 0.63-0.90 
130 345 0.53 0.41-0.67 



By tumor stage (P_^g, = 1 x 10^) 

Carcinoma m s/fu 826 1591 0.67 0.60-0.76 
Invasive carcinoma 134 134 1.38 0.98-1.95 



By age {P_,^,= 0.57) 



<36 
>36 



477 

483 



851 0.71 0.60-0.84 
874 0.76 0.64-0.89 



<4> 



MICA-A5.1 

p= 1.6 X 10 ' Overall 

Heterozygotes 
Homozygotes 



Ca Co OR 95%CI 

949 1689 1.34 1.20-1.50 

454 811 1.45 1.17-1.81 

342 480 1.84 1.46-2.32 



By tumor stage (P_^j, = 0.62) 

Carcinoma m srtu 817 1558 1.33 1.18-1.50 

Invasive carcinoma 132 131 1.46 1.03-2.07 



By age (P.^„ = 0.24) 



<36 
>36 



0.5 1.0 1.5 2.0 

OR (95%CI) 



473 
476 



831 1.44 1.23-1.69 
858 1.26 1.08-1.48 



•<!>P= 3.8x10"' 



1.0 1.21.41.6 2.0 2.4 
OR (95%CI) 



MICA-A4 


Ca 


Co 


OR 


95%CI 


Overall 


949 


1689 


0.80 


0.68-0.94P 


Heterozygotes 


223 


474 


0.79 


0.65-0.95 


Homozygotes 


19 


46 


0.70 


0.41-1.21 


By tumor stage (P. 


het ~ 


0.40) 






Carcinoma in situ 


817 


1558 


0.82 


0.69-0.97 


Invasive carcinoma 


132 


131 


0.65 


0.39-1.10 


By age (P_^^,= 0.76) 








<36 


473 


831 


0.82 


0.65-1.03 


>36 


476 


858 


0.78 


0.62-0.98 



0.4 



MICA-A5 

Overall 

Heterozygotes 
Homozygotes 



0.6 0.8 1.01.2 
OR (95%CI) 



Ca Co OR 95%CI 

949 1689 0.60 0.50-0.72 

163 428 0.60 0.49-0.73 

8 36 0.34 0.16-0.75 



By age (P.^^, = 0.53) 

<36 473 
>36 476 



io-«<i> 



By tumor stage (P.^^, = 0.47) 

Carcinoma in s/Ju 817 1558 0.58 0.48-0.71 
Invasive carcinoma 132 131 0.71 0.42-1.18 



831 0.56 0.44-0.73 
858 0.63 0.48-0.81 



0.2 0.5 1.0 
OR (95%CI) 



Figure 1. Stratified analysis of association between rs9272143, MICA-A5.1, -A4, and -AS and risk of cervical cancer. Unless specified, the odds 
ratios (ORs) and 95% confidence intervals (CIs) of per variant allele (log-additive model) and per genotype were calculated using unconditional 
logistic regression with adjustment of age at recruitment and study design (carcinoma in situ vs. invasive carcinoma) when appropriate. P for 
heterogeneity (P-het) was derived from the Cochran Q test. Squares represent odds ratios; size of the square represents inverse of the variance of 
the log odds ratio; horizontal lines represent 95% confidence intervals; diamonds represent overall estimate; solid vertical lines represent an odds 
ratio of 1; dashed vertical lines represent the overall odds ratios. Ca, case subject; Co, control subject. 



between HLA-DPB2 variant rs3117027 and risk of cervical 
cancer was not replicated in this study. 

SNP rs9272143 is located 4.38 kb upstream of HLA- 
DQAl and 43.19 kb upstream of HLA-DRBl. It has 
recently been identified as a ds-expression quantitative 
trait locus (c!S-eQTL), with the T allele being associated 
with decreased expression of HLA-DRBl as compared to 
the C allele [9]. HLA-DRBl belongs to the HLA class II 
/]-chain paralogs, which encodes the /?-chain of the pep- 
tide-antigen receptor HLA-DR. It is expressed in Lan- 
gerhans cells (LC), the antigen presenting cells of 
squamous epithelia in the cervix, and plays a central 
role in the cell-mediated immune response by presenting 
processed foreign antigens to CD4-I- helper T-lympho- 
cytes [10]. CD4-I- T-cell activation results in the secre- 
tion of a variety of small proteins, or cytokines. Our 



study points to the importance of expression level of 
HLA-DRBl in cervical carcinoma. Impaired class II gene 
expression [11-13] and a reduced number of LC have 
been reported in genital HPV infections [14, 15] and in 
lesions due to HPV [16]. The increased incidence and 
progression of HPV infections in immunosuppressed 
individuals illustrates the critical importance of the 
CD4-I- T-cell-regulated cell-mediated immune response 
in the resolution and control of HPV infection [10, 17]. 
Regression of anogenital warts is accompanied histologi- 
cally by a CD4-I- T-cell-dominated Thl response. Failure 
to develop effective cell-mediated immune response to 
clear or control infection results in a persistent infection 
and, in the case of the oncogenic HPVs, an increased 
probability of progression to CIS and invasive carci- 
noma [17]. Therefore, HLA-DRBl may be involved in 
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the key event of establishing a robust defense against 
HPV infection through antigen presentation to CD4+ 
helper T-lymphocytes. It is, hence, biologically plausible 
that carriers of the C allele of rs9272143, which have 
higher expression level of HLA-DRBl are less susceptible 
to cervical cancer. However, it is unknown whether 
rs9272143 is the pathogenic variant or other functional 
variant(s) which affect(s) the expression of HLA-DRBl 
could be responsible for this signal. Functional studies 
would help to provide insight into it. 

MICA encodes a membrane-bound protein, which acts 
as a ligand to stimulate an activating receptor, NKG2D, 
expressed on the surface of essentially all human natural 
kiUer (NK), y5 T, and CD8+ T cells [18-20]. Nor- 
mally, MICA is constitutively expressed in low levels on 
epithelial cells in the gut and thymus, endothelial cells, 
fibroblasts, and monocytes [21-23]. But it is upregulated 
or expressed de novo in stressed conditions, such as 
during viral and bacterial infections [20, 24, 25], heat 
shock [22], DNA damage response [26], oncogenic 
transformation [18, 19], and in autoimmune conditions 
[27]. MICA serves as signal of cellular stress, and 
engagement of NKG2D by MICA triggers NK cells, and 
costimulates some y5 T cells and antigen-specific CD8^ 
a/? T cells, resulting in a range of immune effector func- 
tions, such as cytotoxicity and cytokine production [21, 
28]. The recognition of the MICA molecule by the 
NKG2D receptor enables immune cells to identify and 
attack infected or transformed cells without the need of 
MHC class I expression or antigen recognition [29]. 
Thus, the MICA/NKG2D interaction is an effective 
mechanism for immunosurveillance. This strong selection 
pressure seems to have led tumor cells to evolve mecha- 
nisms to minimize or avoid the response mediated by 
NKG2D by shedding MICA from the cell surface [30, 
31]. The shedding of MICA has been reported to be 
mediated by metaUoproteinases through proteolytic 
cleavage of the extracellular parts and palmitoylation of 
two cysteine residues in the cytoplasmic tail of MICA 
was found to be necessary for efficient cleavage [31, 32]. 
The shedding of soluble MICA by human tumors not 
only hinders recognition of the MICA-expression tumor 
cells, but also results in systemic downregulation of 
NKG2D on NK and CD8+ T cells, and evasion of 
NKG2D-mediated immune recognition [30, 31]. 

The SNP rs2516448 is located 7.32 kb downstream of 
the MICA gene and its effect is independent of previously 
known associations with HLA alleles/haplotypes [5]. The 
T allele of rs25 16448 that increased the susceptibility to 
cervical cancer is in perfect LD with the MICA-A5. 1 allele, 
which encodes a truncated protein lacking part of the 
TMD and the whole cytoplasmic tail and is most com- 
monly seen in the MICA*008 allele [6]. Suemizu et al. 
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[33] found that the cytoplasmic tail-deleted MICA-A5.1 
gene product was aberrantly transported to the apical sur- 
face of human intestinal epithelial cells instead of the ba- 
solateral surface where the interaction with intraepithelial 
T and NK lymphocytes takes place. Thus, MICA-A5.1 car- 
riers may have an aberrant immunological surveillance by 
NK and T cells. Meanwhile, in contrast to other MICA 
alleles that are shed as truncated soluble species after pro- 
teolysis by metaUoproteinases, the protein translated from 
the MICA-A5.1 allele is released from cells as a mem- 
brane-anchored full-length molecule in exosomes due to 
the lack of the two cysteines required for proteolytic 
shedding. Incubation of NK cells with the MICA-A5.1 
(MICA*008) containing supernatant triggers significantly 
more NKG2D downregulation than the MICA*019 cul- 
ture supernatant. Strikingly, incubation with exosomes 
containing MICA-A5.1 (MICA*008) also impairs NK-cell 
cytotoxicity [34]. Taken together, the preceding evidence 
might explain the result in this study that individuals car- 
rying the MICA-A5.1 allele have a predisposition for their 
infected or transformed cells to escape from attack by 
immune cells. The MICA-A5.1 allele has also been associ- 
ated with different autoimmune diseases and other tumor 
forms [35-42], supporting its role in immune response 
and tumor development. 

In contrast to MICA-A5.1, alleles encoding only four 
(M/CA-A4) or five (MICA-AS) Ala residues in the TMD 
of MICA protein, respectively, were found to confer pro- 
tection against cervical cancer. In particular, a very strong 
protective effect was observed for the A5/A5 genotype. 
Further studies with larger sample sizes are warranted to 
verify this result, given the limited number of MICA-A5 
homozygotes in this study. Interestingly, in breast cancer, 
MICA-A4 and MICA-A5 have also shown a protective 
effect [41, 42]. The mechanism by which MICA-A4 and 
MICA-A5 protect against cervical cancer is yet unknown. 
The nearby HLA-B*0702 has been reported to be associ- 
ated with cervical cancer risk [5, 43, 44]. However, in a 
recent study, we found that the association with HLA- 
8*0702 was actually driven by the joint effects of both 
rs9272143 and MICA-A5.1 (unpublished data), suggesting 
that HLA-B*0702 is unlikely to be the causal variant 
responsible for the association with MICA-A4 or MICA- 
AS. These short tandem repeats are not located in any of 
the extracellular domains and do not directly affect the 
putative binding site of MICA with NKG2D. As the 
amino acids encoded by the microsatellite are situated in 
the TMD of the molecule, it is possible that certain alleles 
provide a more stable anchoring of MICA to the cell sur- 
face and therefore permit better binding to NKG2D and, 
as a consequence, a more efficient NKG2D-mediated 
immunosurveillance [41, 45]. It is also worth noting that 
MICA-A4 is in high LD with the amino acid substitution 
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of Glycine (Gly) by Tryptophan (Trp) at position 14 in 
the al domain and MICA-A5 is in high LD with the 
amino acid substitution of Gly by Serine (Ser) at position 
175 in the a2 domain of MICA in the IMGT/HLA 
database [46]. The top surface of the MICA al-oi2 plat- 
form has been found to interact directly with NKG2D 
[47, 48]. Further studies are warranted to determine 
whether these amino acid changes could affect the bind- 
ing affinity of MICA to the NKG2D receptor and whether 
they are responsible for the associations with MICA-A4 
and MICA-A5. 

The SNP rs3 117017 resides in the pseudogene HLA- 
DPB2, but close to the functional gene HLA-DPBl, which 
encodes the /?-chain of the peptide-antigen receptor HLA- 
DP. In contrast to our previous GWAS of cervical cancer, 
we did not observe a significant association between this 
SNP and risk of cervical cancer. It is noteworthy that the 
effect size of rs3117027 variant was much smaller 
(OR = 1.25) than the other two hits and the P value was 
close to genome-wide significance threshold in the initial 
GWAS discovery cohort (P = 3.1 x 10^*). It is possible 
that this study lacks sufficient statistical power to detect 
its modest effect on susceptibility for cervical cancer. Fur- 
ther studies with larger sample sizes are needed to draw a 
firm conclusion. On the other hand, the subjects enrolled 
in this study are from northern Sweden only, while the 
subjects in the GWAS were a national collection. The lack 
of replication may indicate no direct association between 
rs3117027 and cervical cancer. One cannot dismiss the 
possibility, however, that the LD between rs3117027 and 
the putative causal variant varies between ethnically dis- 
tinct populations. 

This study has several limitations. First, the number of 
invasive cervical cancer cases (134) included in the study 
is modest. Although we have observed evidence for heter- 
ogeneity by tumor stage for rs9272143, the statistical 
power is insufficient to draw any firm conclusion. Future 
studies with larger numbers of invasive cancers are war- 
ranted to validate this finding. Second, except for age, 
information on other risk factors for cervical cancer such 
as HPV infection, parity, oral contraceptive use, and 
tobacco smoking [49], which might modify the effects of 
the susceptibility loci, was not available in our study. Pos- 
sible interactions between the susceptibility loci and these 
risk factors should be thoroughly investigated in future 
studies. 

In summary, associations identified in our previous 
GWAS with rs9272143 in the MHC class II region, as well 
as rs2516448 and the MICA-A5.1 allele in the class I 
region were replicated in a northern Swedish population, 
providing credible evidence that these genetic variants 
influence susceptibility to cervical cancer. We also identi- 
fied a reduction in risk associated with the MICA-A4 and 
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MICA-A5 alleles. Our results do not support previous 
suggestion that HLA-DPB2 variant rs3 117027 is positively 
associated with cervical cancer. The association with 
rs2516448 seems to be driven by the functional allele 
MICA-A5.1. However, it is unknown whether other allele 
combinations not measured by this study could be 
responsible for the signals of rs9272143 as well as MICA- 
A4 and MICA-A5 alleles. Further functional studies are 
warranted to identify the causal variant (s) responsible for 
these signals and their functional effect (s). 
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